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ABSTRACT 



The Fe unresolved transition array (UTAs) produce prominent features in the 
~ 15 — 17 A wavelength range in the spectra of Active Galactic Nuclei (AGN). 
Here we present new calculations of the energies and oscillator strengths of inner- 
shell lines from Fexiv, Fexv, and Fexvi. These are crucial ions since they are 
dominant at inflection points in the gas thermal stability curve, and UTA exci- 
tation followed by autoionization is an important ionization mechanism for these 
species. We incorporate these, and data reported in previous papers, into the 
plasma simulation code Cloudy. This updated physics is subsequently employed 
to reconsider the thermally stable phases in absorbing media in Active Galactic 
Nuclei. We show how the absorption profile of the Fexiv UTA depends on den- 
sity, due to the changing populations of levels within the ground configuration. 

Subject headings: atomic data — atomic processes — line: formation — X-rays: 
galaxies — galaxies: active — methods: numerical 
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Introduction 



X-ray observations of active galactic nuclei (AGN) have revealed the existence of a 
broad absorption feature in the 15-17 A wavelength range. This feature arises from an 
unresolved transition array (UTA) due primarily to 2p — > 3d inner-shell absorption lines in 



iron ions with an open M shell (Fel - Fexvi), and was first identified by Sako et al. (2001) 
in the XMM-Newton spectrum of the quasar IRAS 13349+2438. Since then, the UTA has 



been detected in numerous AGN (see, for example, Longinotti et al. (2010), Lestinsky et al. 



(2009) and references therein). It is believed to arise in the warm, relatively low ionization, 



absorbing gas surrounding the central supermassive black hole in the AGN. 



Behar et al. (2001) first noted the diagnostic potential of the UTA, and showed that 



their central wavelengths, absorption profiles and equivalent widths can provide information 
on the ionization structure, kinematics and column densities of the warm absorbing material 
in the AGN. These authors also produced atomic data required for the modeling of the 
UTA, including energy levels (wavelengths) and oscillator strengths for the 16 iron charge 
states Fel through Fexvi, calculated with the multiconfiguration, relativistic hullac 



computer code (Bar-Shalom et al. 2001). Since then, several authors have calculated 



wavelengths and oscillator strengths for the UTA. These include Kisielius et al. (2003), who 



produced atomic data for Fexv and Fexvi using the configuration interaction code Civ3 



(Hibbert 1975) with the inclusion of relativistic effects by adding Breit-Pauli operators 



to the Hamiltonian (Hibbert et al. 1991). More recently, Gu et al. (2006) calculated 



results for FeVl through Fexvi employing second-order many-body perturbation theory 



(Lingren 1974). Furthermore, Beiersdorfer et al. (2012) used the relativistic multi-reference 



M0ller-Plesset perturbation theory to calculate the energy levels of Fexvi, including 
those of the autoionizing levels with a hole state in the L shell. They obtained very good 
agreement between the calculated L-shell transition wavelengths and those from recent 



-4- 



laboratory measurements. 

Here we present extended calculations of UTA energies and radiative transition rates 
for Fexiv, Fexv and Fexvi. These are combined with other improvements in the plasma 



simulation code Cloudy (Ferland et al. 2013) to generate a range of models of absorbing 



clouds in AGN. Calculations were performed with r6433 of the Cloudy development branch, 
and will be part of the 2013 release of the code. We also perform a thermal stability analysis 
to identify the thermally stable phases where clouds can exist, and show that Fexiv, Fexv 
and Fexvi probe the high ionization end of the warm stable phase. 



2. New atomic data for Fexiv, Fexv and Fexvi transitions 

Previously, a comparatively small set of inner-shell photoexcitation data for Fe XV and 



Fe xvi was reported by Kisielius et al. (2003). In that work energy levels were identified 



and electric dipole transition wavelengths, oscillator strengths and transition probabilities 
determined. Specifically, transitions involving photoexcitation from the inner 2p shell to 
the outer n = 3 shell were investigated using the relativistic Breit-Pauli approach, and an 
assessment of the accuracy of the derived atomic data was given. In the present work, 
the data set has been significantly extended. First of all, we add line data involving 
the transitions from the 2s shell since their wavelengths lie within the same UTA range. 
Furthermore, we determine excitation (both from 2s and 2p shells) to the n = 3,4, 5 shells. 
Finally, photoexcitation data for the same type of transitions in Fe xiv ions are derived. 

In our current work we use the configuration interaction (CI) approximation for atomic 



calculations. Particularly, the well-established computer code Civ3 of (Hibbert 1975) is 
exploited to determine multiconfiguration eigenfunctions and eigenvalues of considered 
energy levels. After this initial step, we derive atomic data such as electric dipole (El) 
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radiative transition wavelengths A, oscillator strengths / and transition probabilities (rates) 
A for the transitions from inner n = 2 shell of Fe xiv, Fe XV and Fe xvi ions. The adopted 
CI method is suitable for production of very accurate data, where the correlation effects 
play a vital role in theoretical data accuracy. 



Data are generated in the relativistic Breit-Pauli approach (Hibbert et al. 1991) for all 
lines arising from 2s and 2p shell electron transitions to the n = 3, 4, 5 shells. Combining 
these two approaches, we can reliably determine LS'J-energy levels and, consequently, the 
El transition data for fine-structure levels. When CI wavefunctions for the fine-structure 
levels of the ground and excited states are determined, the electric dipole transition 
absorption oscillator strength / is given by 



f = 2/3g~ 1 AE ij S 2 (1) 

where g>j = 2Jj + 1 is the statistical weight of the initial (ground) level, AEij is the energy 
difference between the upper and lower level, and S is transition line strength (matrix 
element). The corresponding electric dipole line emission transition probability A (in s _1 ) 
is given by 

A = 2.142 x lO 10 ^- 1 (AEij) 3 S 2 (2) 
where gj = 2Jj + 1 is the statistical weight of the excited-state level j. 

A more detailed description of the approximations used in the present work may be 



found in Kisielius et al. (2003), where inner 2p shell photoexcitation was considered for 



Fexv and Fexvi. The approach applied in the current paper for a much wider set of lines 



and Fe ions is very similar. Kisielius et al. (2003) have considered in detail the accuracy 



of determined results by applying mainly the internal indicators such as convergence of 
transition wavelengths with increasing the size of CI wavefunction expansion or agreement 
between length and velocity forms of oscillator strength values. They have estimated that 
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2p — 31 transition wavelength accuracy was 0.2 % whereas oscillator strength accuracy was 
assessed as 10% for the lines having / > 5 x 10~ 5 . The conclusions on the accuracy of the 



atomic data made by Kisielius et al. (2003) are also applicable to the present results. 



2.1. Fexvi lines 

The ground level of the Fe 15+ ion is ls 2 2s 2 2p 6 3s 2 Sy 2 . Its CI wavefunction expansion 
was chosen to include no more than two electrons excited from the 2p and 3s shells in the 
case of the lines arising from the 2p shell, and no more than two electrons excited from 2s 
and 3s shells for the lines arising from the 2s shell. The significant difference between the 
calculation of n = 3 and 4 lines lies in the fact that we use correlation radial orbitals 4s, 4p, 
4d, 4/ in the former case to improve the calculated wavelength accuracy (for more details 
see Kisielius et al. 2003). In the latter case these orbitals represent real states while the 



configuration sets remain the same. 



2.1.1. Excitation from the 2p shell 

Following this approach, the ground state CI wavefunction expansion for the transition 
arrays 2p — 31 and 2p — Al was the same and included 2p 6 3l,2p 6 4l (I = 0,2), 2p 5 3s3p, 
2p 5 3p3d, 2p 5 3l4l', 2p 5 4l4l' (I = 0, 1, 2; /' = 0, 1, 2, 3) even-parity configurations, coupled 
to the 2 S, 2 P, 4 P, A D LS-terms. This CI expansion generates 94 LSJ-\eve\s with total 
orbital quantum number J = 1/2. Due to selection rules for El transitions, the excited 
configuration with the vacancy in the 2p shell can only have levels with J = 1/2 and 
J = 3/2. We use a CI wavefunction expansion consisting of odd-parity configurations 2p 6 3p, 
2p 6 4p, 2p 5 3l 2 (I = 0,1,2), 2p 5 3s3d, 2p 5 3M' (I = 0,1,2;/' = 0,1,2,3), 2p 5 4l 2 (l = 0,1,2,3), 
2p 5 4l4l' (I = 0, = 2,3) coupled to the 2 S, 2 P, 4 P, 4 D LS-terms for J = 1/2 and to 
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the 4 S, 2 P, 4 P, 2 D, 4 D, 4 F terms for J = 3/2. Consequently, we consider 93 levels with 
J = 1/2 and 147 levels with J = 3/2 for the excited configurations having electrons in the 
outer n = 3 and 4 shells and the vacancy in the inner 2p shell. 

Very similar configuration sets were used to obtain atomic data for the transitions 
from the 2p shell to the outer 5/ shell. The ground state wavefunction expansion included 
2p 6 nZ (n = 3,4,5;/ = 0,2), 2p 5 3s3p, 2p 5 3p3d, 2p 5 3/5/' (/ = 0,1,2;/' = 0,1,2,3), 2p 5 5/5/' 
(/ = 0,1,2,3;/' = 0,1,2,3) even-parity configurations coupled to the 2 S, 2 P, 4 P, 4 D 
LS'-terms. This yields 95 LSJ-levels having J = 1/2. 

The same (as in n = 3,4 case) terms for the excited configuration levels with 
J = 1/2 and J = 3/2 were considered, while the CI expansion included configurations 
2p 6 np (n = 3,4,5), 2p 5 3/ 2 (/ = 0,1,2), 2p 5 3s3d, 2p 5 3/5/' (/ = 0,1,2;/' = 0,1,2,3) 2p 5 5/ 2 
(/ = 0, 1, 2, 3), 2p 5 5s5d, 2p 5 5p5f. For J = 1/2, the CI expansion included 94 levels whereas 
for J = 3/2 the number was 148. Levels arising from the configurations with 5g electrons 
were not included in the CI expansion because such lines are very weak compared to other 
n = 5 transitions. 



2.1.2. Excitation from the 2s shell 

A similar approach for the configuration sets to that adopted for the 2p shell was 
taken when calculating atomic data for the lines corresponding to the transitions from the 
inner 2s shell to the outer n = 3, 4, 5 shells. For the n = 3, 4 lines, the ground state CI 
expansion included 56 levels with J = 1/2 from the configurations 2p 6 nl (n — 3, 4; / = 0, 2) 
2p 5 3s3p, 2p 5 3p3d, 2s2p 6 3l 2 (/ = 0,1,2), 2s2p 6 3/4/' (/ = 0,1,2;/' = 0,1,2,3), 2s2pHl 2 , 
(/ = 0, 1,2,3), 2s2p 6 4s4d, 2s2p 6 4p4f coupled to the 2 S, 2 P, 4 P, 4 P> terms. The upper-state 
configuration set consisted of 2p 6 3p, 2p 6 Ap, 2p 6 Af, 2s2p 6 3s3p, 2s2p 6 3p3d, 2s2p 6 nlAl' 
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(n = 3, 4; / = 0, 1, 2; V = 0, 1, 2, 3) configurations coupled to the terms 2 S, 2 P, 4 P, 4 P>, 6 D, 
6 P for the J = 1/2 levels, and to the 2 P, 2 P>, 4 S, 4 P, 4 P>, 4 P, 6 P, 6 P>, 6 P, 6 G terms for 
J = 3/2. This selection produced 38 and 58 levels for J = 1/2 and 3/2, respectively. 

For the lines corresponding to the transitions from the inner 2s shell to the outer 
n = 5 shell, the ground state wavefunction CI expansion included configurations 2p 6 nl 
(n = 3,4,5;/ = 0,2), 2p 5 3s3p, 2p 5 3p3d, 2s2p 6 3l 2 (/ = 0,1,2), 2s2p 6 3s3d, 2s2p 6 3/5/' 
(I = 0,1,2;/' = 0,1,2,3), 2s2p 6 5/ 2 (/ = 0,1,2), 2s2p 6 5s5d, 2s2p 6 5p5f coupled to the 
terms 2 S, 2 P, 4 P, 4 D. This generates a total of 57 levels with J = 1/2. For the upper 
state configuration set we included configurations 2p 6 nl (n = 3,4,5;/ = 1,3), 2s2p 6 3s3p, 
2s2p 6 3p3d, 2s2p 6 n/5/', (n = 3, 4; I = 0, 1, 2; V = 0, 1, 2, 3) which were coupled to the terms 
2 S, 2 P, 4 P, 4 P, 6 P, 6 P, producing 39 levels with J = 1/2. For J = 3/2 we considered the 
same terms as in the n = 4 lines case, producing 59 levels. 



2.1.3. Transition data 

We employ multiconfiguration wavefunctions determined after the Hamiltonian 
matrix diagonalization to calculate transition wavelengths A, oscillator strengths /, and 
transition rates A for excitation from the inner n = 2 shell to the outer n = 3, 4, 5 shells. 
Calculated line data cover the wavelength region A = 12.606 — 14.096 A for 2s — 31 
transitions (a total of 17 lines), A = 10.145 — 11.221 A for 2s — 4/ transitions (46 lines) and 
A = 9.545 — 10.333 A for 2s — 5/ transitions (46 lines). Considering transitions from the 2p 
shell, the wavelength range was A = 13.594 — 17.246 A for the excitation of the n = 3 shell 
(45 lines), A = 11.316 - 12.934 A for the n = 4 shell (111 lines), and A = 10.525 - 11.689 A 
for the n = 5 shell (111 lines). 

In Table[T]we present data for Fexvi transitions from the 2s and 2p shells to n = 3, 4, 5. 
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Table 1. Transitions 2p — rijd and 2s — rijp from the L shell to rij = 3, 4, 5 shells in Fe 15+ . 
Column riili denotes the initial shell, rij the final shell for optical electron transition, and 
gj the statistical weights for the initial and final levels, A the line wavelength (in A) and 

the absorption oscillator strength. 
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Note. 




Table 


[l]is 


published in its 



entirety in the electronic edition of the 
Astrophysical Journal. Only a portion 
is shown here for guidance regarding its 
form and content. 
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Fig. 1. — A comparison of the wavelengths A and the absorption oscillator strengths fij for 



the strongest lines in Fexvi. Dots represent our data, crosses are the results of Gu et al. 



(2006) and stars are those from Behar et al. (2001). The data at A = 15 A represent 2p — 3d 



UTA, at A = 14A represent 2s - 3p UTA, at A = 12.5A represent 2p - Ad UTA and at 
A = 11. 5A represent 2p — 5d UTA. The Behar et al. (2001) data are given for an /-value 



averaged wavelength A av and the total oscillator strength /y of the line group. See the 
electronic edition of the Journal for a color version of this figure. 

We include only transitions with / > 0.1, whereas the more complete online version of the 
table contains lines having / > 0.0001. Weaker lines are not presented since they do not 
affect the ionization balance nor the final spectra. More details on data trimming are given 
in Sect. 13.21 



In Fig. [T] we compare the strongest lines from our calculation with data derived using a 
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many-body perturbation theory (MBPT, Gu et al. 2006), and with the Behar et al. (2001) 



results obtained using the multiconfiguration relativistic HULLAC computer package. The 
Behar et al.] ( |2001| ) data are given for an /-value averaged wavelength A av and the total 
oscillator strength ^ fij of the line group, and therefore there is a sizable difference in the 
number of lines and in the values of /. Note that the point at A = 17.296A represents 
transitions from the 2p to 3s shells which are not considered in our calculations. When 



comparing our results with those of Gu et al. (2006), one can see quite good agreement both 



for wavelengths and oscillator strengths. There is some systematic shift towards shorter 



wavelengths in our data, as discussed in Kisielius et al. (2003). In addition, a significant 



difference with Gu et al. (2006) is the appearance in our results of a group of lines at 
A = 11.5 — 12.5 A which represent transitions to n = 4 and n = 5 states. Some of these 
have large /-values and are important to include in atomic models. 

The primary goal of the calculations presented here is to create an extensive set 
of inner-shell transition data to ensure that the total excitation and ionization rate is 
computed as accurately as possible. In a later section we show how this new data changes 
the ionization distribution of gas near an AGN. The line wavelengths we quote, which come 
from our atomic structure calculation, are no better or worse than the wavelengths quoted 
in previous studies. Indeed, the scatter in the wavelengths represents the uncertainty. 
Extensive sets of laboratory measurements of energy levels would be needed to improve 
the level energies and resulting wavelengths. Although a few experiments have been done 



(Brown et al. (2001) and Simon et al. (2010)), the extensive laboratory data needed to 



significantly improve the line wavelengths do not now exist. 
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2.2. Fexv lines 

The ground level of the Fe 14+ ion is ls 2 2s 2 2p 6 3s 2 1 Sq. Its CI wavefunction expansion 
was chosen to include no more than two electrons virtually excited from the 2p and 3s 
shells in the case of the lines arising from the 2p shell, and no more than two electrons 
excited from the 2s and 3s shells for the lines arising from the 2s shell. We adopt the 
same configuration set both for the 2p — 31 and the 2p — 4/ line calculations, but different 
radial orbital sets are used. For the n = 3 line calculation, we employ the correlation radial 
orbitals 4s, 4p, 4c/, 4/, whereas the corresponding radial orbitals represent real states in the 
n = 4 line calculations. This method contributes to an increased accuracy of the calculated 
wavelengths for n = 3 lines which are much stronger than the n = 4 ones. 



2.2.1. Excitation from the 2p shell 

For excitation from 2p to the n = 3,4 shells, a CI wavefunction expansion was 
made of the 2p 6 3l 2 {I = 0,1,2), 2p 6 3s3d, 2p 5 3s 2 3p, 2p 5 3p 3 , 2p 5 3p3d 2 , 2p 5 3s3p3d 
configurations containing two or three electrons within the n = 3 shell. The even-parity 
configurations 2p 6 3/4/' (I = 0,1,2;/' = 0,1,2,3), 2p 5 3/ 2 4/',(/ = 0,1,2;/' = 1,3), 2p 5 3/3/'4/" 
(/ = 0, 1; /' = 0, 1, 2; /" = 0, 1, 2, 3) with a single electron virtually excited to the n = 4 shell 
were added to that set as well as even-parity configurations 2p 6 Al 2 (/ = 0, 1, 2, 3), 2p6AsAd, 
2pHp4f, 2p 5 3pAl 2 (/ = 0,1,2,3), 2p 5 3/4/'4/" (/ = 0,1,2;/' = 0,1,2;/" = 1,2,3) with two 
electrons excited to the n = 4 shell. In the case of inner 2p shell excitation, we consider the 
1 S e , 3 P e and 5 D e terms which produce the J = fine-structure levels. Construction of such 
a basis produces 408 configuration state functions (CSFs) in the CI wavefunction expansion 
for the lower state of the Fe 14+ ion. 

Due to the electric dipole transition selection rules, we have to consider fine-structure 
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levels with J = 1 in the upper (odd) state of this ion. In the same way as for 
the lower (even) state, we construct three subsets of configurations in wavefunction 
expansion for the upper (odd) state. The CI wavefunction expansion consists of 
the 2p 6 3s3p, 2p 6 3p3d, 2p 5 3s 2 3d, 2p 5 3s3p 2 , 2p 5 3s3d 2 , 2p 5 3p 2 3d, 2p 5 3d 3 odd-parity 
configurations which have two or three electrons in the n = 3 shells. The odd-parity 
configurations 2p 6 3/4/' (/ = 0, 1, 2; /' = 0, 1, 2, 3), 2p 5 3/ 2 4/' (I = 0, 1, 2; V = 0, 2), 2p 5 3/3/'4/" 
(/ = 0,1; V = 0,1,2;/" = 0,1,2,3) with all possible electron distributions having a 
single electron virtually excited to the n = 4 shell and the same parity configurations 
with two electrons in open n = 4 shells 2p 6 4/4/' (I = 0,1,2;/' = 1,2,3), 2p 5 3/4/' 2 
(/ = o, 2; Z' = 0, 1,2,3), 2p 5 3/4/ / 4/" (/ = 0, 1, 2; V = 0, 1, 2; I" = 1, 2, 3 were added to that 
expansion. 

In this construction of the CI wavefunction basis, we restricted our set by allowing no 
more than two electrons in the excited n = 4 shell. Here we consider the 3 S°, 1 P°, 3 P°, 5 P°, 
3 D°, 5 D°, 5 F° terms which can produce fine-structure levels with J = 1. As a result, we 
have the CI wavefunction expansion with 1094 CSFs for the upper states of the Fe 14+ ion. 

A very similar configuration set was applied when considering 2p — 5l transitions. In this 
instance, the configurations with n = 4 electrons were replaced by configurations with n — 5. 
The only significant difference was the addition of the 2p 6 3/4/' (/ = 0, 1,2;/' = 0, 1,2,3) 
even-parity configurations for in CI expansion of the lower state of Fe 14+ , and the 2p 6 3/4/' 
(/ = 0,1,2;/' = 0,1,2,3) odd-parity configurations for the excited state. Configurations 
with 5g electrons were not included in the CI wavefunction expansion because transitions 
to such configurations are very weak. The total number of fine-structure levels with J = 
in this case was 413 for the lower configuration and 1107 for the upper configuration levels 
with J = 1. 
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2.2.2. Excitation from the 2s shell 

While determining the transition array 2s — 3/, 4/ lines, for the CI expansion we have 
adopted the same subset (as in 2p — 31,41 case) of configurations with closed 2s 2 and 
2p 6 electron shells and complemented it with configurations having one electron virtually 
excited from the 2s shell. Specifically for the ground level, the even-parity configurations 
2s2p 6 3s 2 3d, 2s2p 6 3p 2 3d, 2s2p 6 3s3p 2 , 2s2p 6 3s3d 2 , 2s2p 6 3d 3 , 2s2p 6 3/ 2 4/' (/ = 0,1,2;/' = 0,2), 
2s2p 6 3/3/'4/" (/ = 0,1; I' = 1,2;/" = 0,1,2,3), 2s2p 6 3/4/' 2 (/ = 0,2;/' = 0,1,2,3), 
2s2p 6 3/4/'4/" (/,/' = 0, 1,2;/" = 1,2,3), were included in our CI wavefunction expansion. 
These configurations were coupled in the 1 S°, 3 P° and 5 D° non-relativistic LS'-terms and 
produced 165 fine-structure levels with total orbital quantum number J = 0. 

The construction of the configuration set for the upper state is based on the same 
principles. We employ configurations with closed 2s 2 and 2p 6 electron shells from the 
2p — 3/, 4/ array calculation and add the configurations with an open 2s shell. Following this 
method, we include the configurations 2s2p 6 3s 2 3p, 2s2p 6 3p3d 2 , 2s2p 6 3p 3 , 2s2p e 3s3p3d with 
the 2s electron virtually excited to the n = 3 shell, the odd-parity configurations 2s2p 6 3/ 2 4/' 
(/ = 0, 1, 2; V = 1, 3), 2s2p 6 3/3/'4/" (/ = 0, 1, 2; V = 1, 2; I" = 0, 1, 2, 3) with the 2s electron 
virtually excited to the n = 4 shell, and the odd-parity configurations with two electrons in 
the n = 4 shell 2s2p 6 3/4/'4/", (/ = 0, 1, 2; V = 0, 1, 2; /" = 1, 2, 3). Coupled to the 3 S°, 1 P°, 
3 P°, 5 P°, 3 D°, 5 D°, 5 F° terms, these configurations generate total 417 fine-structure levels 
with J = 1. 

For the transition array 2s — 5/, the wavefunction CI expansion was similar to that 
for the 2s — 4/ lines with the n = 4 shell electrons replaced by n = 5. The only addition 
was even-parity configurations 2p 6 3lAl' (/ = 0,1,2;/' = 0,1,2,3 with one electron in the 
4/ shell, included in the lower-state CI wavefunction expansion set, and the odd-parity 
configurations 2p 6 3/4/' (/ = 0, 1, 2; /' = 0, 1, 2, 3 in the upper-state wavefunction expansion. 
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This CI expansion selection produces 170 fine-structure levels with J = for the ground 
state and 430 levels with J = 1 for the excited state. 



2.2.3. Transition data 

Similarly to the Fe xvi case, we obtained multiconfiguration wavefunctions to produce 
photoexcitation data for the Fe XV UTA lines. The calculated line data cover the wavelength 
region A = 11.592-14.199 A for 2s - 31 transitions (a total of 35 lines), A = 9.620-11.423 A 
for 2s - 41 transitions (147 lines) and A = 9.545 — 10.333 A for 2s - 51 transitions (147 lines). 
Considering transitions from the 2p shell, the wavelength range was A = 12.475 — 16.266 A 
for the excitation of the n = 3 shell (88 lines), A = 10.630 - 13.244 A for the n = 4 shell 
(384 lines), and A = 9.998 - 12.017 A for the n = 5 shell (384 lines). 

In Table [2] we present a sample of the line data for Fe XV. Transition wavelengths A 
and oscillator strengths £y for excitation from the 2s and 2p shells to the outer n = 3,4, 5 
shells are presented. Only a small portion of lines with / > 0.4 is included, with the more 
complete online version of the table containing all lines which have / > 0.0001. It is evident 
from Table [2] that the lines representing transitions to the n = 3,4 shells are sufficiently 
strong as to affect spectral formation, and therefore they cannot be excluded from UTA 
line data sets. 

Fig. [2] compares the strongest lines from our calculations with the results derived 



using the MBPT approach (Gu et al. 2006) and with the relativistic data from Behar et al. 



(2001). In general, the agreement between all three sets of data is satisfactory. There is a 



systematic shift towards shorter wavelengths in our data set as noted for the Fe xvi lines. 
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Table 2. Transitions 2p — rijd and 2s — njp from the L shell to the rij = 3, 4, 5 shells in 
Fe 14+ . Column riiU denotes the initial shell, rij the final shell for optical electron transition, 
Qi and gj the statistical weights for initial and final levels, A the line wavelength (in A) and 

fij the absorption oscillator strength. 
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Fig. 2. — A comparison of the wavelengths A and absorption oscillator strengths for the 



strongest lines in Fexv. Dots represent our data, crosses the results of Gu et al. (2006) 



and stars those from Behar et al. (2001). The data at A = 15. 5A represent 2p — 3d UTA, 



those at A = 14. lA represent 2s — 3p UTA, and those at A = 12. 7A represent 2p — 4d UTA. 



The Behar et al. (2001) data are given for an /-value averaged wavelength A av and the total 



oscillator strength ^ of the line group. See the electronic edition of the Journal for a 
color version of this figure. 

2.3. Fexiv lines 



The ground configuration of the Fe 13+ ion is ls 2 2s 2 2p 6 3s 2 3p. It has two fine-structure 
levels, 2 P°/ 2 an d 2 P^/2i separated by 2.33741 eV. Consequently, while considering 
electric-dipole transitions, we must determine energy levels with total angular momenta 
J — 1/2, 3/2, 5/2 of the excited configuration of the Fe 13+ ion. 
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We apply real radial orbitals for Is, 2s, 2p, 3s, 3p, 3d electrons and correlation ones 
for the virtually excited 4s, 4p, Ad, 4/ electrons. This is consistent with the methods 



used to determine transition data for Fexvi and Fexv (see Kisielius et al. (2003)), for 
photoexcitation from the 2p and 2s shells of the ground configuration to the levels of the 
excited n = 3 configurations. Furthermore, for the excitation of the n = 4 and n = 5 lines, 
we adopt only real radial orbitals to describe 4/ and 5/ (I = 0, 1, 2, 3) electrons. 



2.3.1. Excitation from the 2p shell 

We employ the same configuration set both for calculations of 2p — 31 and 2p — 4/ 
transition data. The CI expansion for the lower state includes configurations which have 
one or two electrons virtually excited from the n = 3 shell. These are complemented 
with configurations which have the 2p electron moved to the n = 3 shell. The CI 
wavefunction expansion for the ground state was composed of the configurations with 
three or four electrons in n = 3 shell 2p 6 3s 2 3p, 2p 6 3p3d 2 , 2p 6 3p 3 , 2p 6 3s3p3d, 2p 5 3s 2 3p 2 , 
2p 5 3p 2 3d 2 , 2p 5 3s3d 3 , 2p 5 3p 4 , 2p 5 3d 4 , 2p 5 3s3p 2 3d, the odd-parity configurations with 
one electron in virtually-excited n = 4 shell 2p 6 3l 2 Al' (I = 0, 1, 2; V = 1, 3) 2p 6 3l3l'4l" 
I — 0, 1; I' — 1, 2; I" = 0, 1, 2, 3), and the odd-parity configurations with two electrons in 
n = 4 shell 2p 6 3/4/'4/ ,/ 1 = 0,1, 2; I' = 0, 1, 2; I" = 1, 2, 3), 2p 6 3p4/ 2 (I = 0, 1, 2, 3). In the case 
of the J = 1/2 levels, we consider the 2 5*, 2 P, 4 P, 4 D, e D, 6 F odd-parity LS-terms making 
up 336 CSFs in the CI wavefunction expansion, while for the J = 3/2 fine-structure levels, 
the odd-parity LS-terms 2 P, 2 D, 4 S, 4 P, 4 D, 4 F, 6 P, 6 D, 6 F, 6 G produce 553 CSFs. 

For the upper state levels, the CI wavefunction expansion for the 2p — 4/ lines 
calculation is constructed from even-parity configuration sets. The configurations with 
electrons in n = 3 shell 2p 6 3s 2 3d, 2p 6 3s3p 2 , 2p 6 3s3d 2 , 2p 6 3p 2 3d, 2p 6 3d 3 the even- 
parity configurations with one electron in n = 4 shell 2p 6 3/ 2 4/' (/ = 0,1,2;/' = 0,2), 
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2p 6 3/3/'4/", (/ = 0,1; V = 1,2; I" = 0,1,2,3) the even-parity configurations with two 
electrons virtually excited into the n = 4 shell, namely 2p 6 3/4/' 2 (/ = 0,2;/' = 0,1,2,3), 
2p 6 3/4/'4/" (I = 0,1,2;/' = 0,1,2;/" = 1,2,3) and the even-party configurations with 
a vacancy in the inner 2p shell: 2p 5 3s 2 3p3d, 2p 5 3s3p3d 2 , 2p 5 3s3p 3 , 2p 5 3p3d 3 , 2p 5 3p 3 3d, 
2p 5 3s 2 3/4/' (/ = 1,2;/' = 0,1,2,3), 2p 5 3s3/ ,2 4/" (/ = 0,1;/' = 1,2;/" = 0,1,2,3), 2p 5 3/ 3 4/' 
(/ = 1,2;/' = 0,1,2,3), 2p 5 3/ 2 4/'4/", (/ = 0,1,2;/' = 0,1,2;/" = 1,2,3), 2p 5 3/3/'4/" 2 
(I = 0, 1; V = 1, 2; /" = 0, 1, 2, 3), 2p 5 3/3/'4/"4/"' (/ = 0, 1; V = 1, 2; /" = 0, 1, 2; /"' = 1, 2, 3). 

These configurations are bound to the 2 S, 2 P, 4 P, 4 L>, 6 P, 6 L>, 6 P LS-terms for 
J = 1/2 and produce 4099 CSFs in the CI wavefunction expansion. The terms 2 P, 2 D, 4 S, 
4 P, 4 L>, 4 P, 6 P, 6 L>, 6 F, 6 G generate a total of 7049 fine-structure levels for J = 3/2. For 
the total angular momentum J = 5/2, the above configurations are bound to the 2 D, 2 F, 
4 P, 4 P>, 4 P, 4 G, 6 S, 6 P, 6 P, 6 P, 6 G, 6 H non-relativistic terms and produce 8211 levels. 

The CI expansion for the lower state in the case of 2p — 5/ consists of configuration 
sets similar to those used in 2p — 31, 4/ lines calculation. However, the four configuration 
complexes with 4/ (/ = 0,1,2,3) electrons are replaced by the configurations with 5/ 
electrons. Furthermore, the CI wavefunction expansion is extended by an additional 
set of the odd-parity configurations 2p 6 3/ 2 4/' (/ = 0,1,2;/' = 1,3) and 2p 6 3/3/'4/" 
(/ = 0, 1; /' = 1, 2; I" = 0, 1, 2, 3) with one electron in an outer n = 4 shell. The same (as 
in case of transitions to n — 3, 4) non-relativistic LS-terms produce 359 CSFs in the CI 
wavefunction expansion for J = 1/2, and 591 CSFs for J = 3/2. 

Similarly, the configuration sets with 4/ electrons are replaced by the configuration 
sets with 5/ electrons for the upper states. This set is extended with the even-parity 
configurations having one electron in the n = 4 shell: 2p 6 3/ 2 4/' (/ = 0, 1,2; I' = 0,2) and 
2p 6 3/3/'4/" (/ = 0, 1; V = 1, 2; /" = 0, 1, 2, 3). As in the previous case, even-parity LS-terms 
produce 4326, 7421 and 8612 CSFs in the CI wavefunction expansion for total angular 
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momentum J — 1/2, 3/2 and 5/2, respectively. 

2.3.2. Excitation from the 2s shell 

For the lines representing transitions from the inner 2s shell to the valence n = 3 or 
n = 4 shells, the CI wavefunction expansion configurations for the lower state consists of 
the same configuration complexes for both cases. The first are configurations with outer 
electrons in n = 3 shell 2p 6 3s 2 3p, 32p 6 p 3 , 2p 6 3p3d 2 , 2p 6 3s3p3d and those with 2s vacancy 
2s2p 6 3s 2 3p3d, 2s2p 6 3s3p3d 2 , 2s2p 6 3s3p 3 , 2s2p 6 3p 3 3d, 2s2p 6 3p3d 3 . These configurations 
are complemented with odd-parity configurations with one and two electrons in the 
4/ shell: 2p 6 3/ 2 4/' (I = 0,1,2;/' = 1,3), 2p 6 3/3/'4/" (I = 0,1;/' = 1,2;/" = 0,1,2,3), 
2p 6 3/4/'4/" (/ = 0,1,2;/' = 0,1,2;/" = 1,2,3), 2p 6 3p4l 2 (/ = 0,1,2,3), 2s2p 6 3/ 2 3/'4/" 
(/ = 0, 1; /' = 1, 2; I" = 0, 1, 2, 3), 2s2p 6 3/3/ ,2 4/" (/ = 0, 1; /' = 1, 2; /" = 0, 1, 2, 3), 2s2j9 6 3/ 3 4/' 
(/ = 1, 2; /' = 0, 1, 2, 3), and 2s2p 6 3s3p3dAl (/ = 0, 2). 

In the case of the J = 1/2 levels, the above configurations are bound to the 2 S, 2 P, 4 P, 
4 _D, 6 D, 6 F odd-parity LS'-terms which yield a total of 660 CSFs in the CI wavefunction 
expansion. For the J = 3/2 levels represented by the 2 P, 2 D, 4 S, 4 P, 4 D, 4 F, 6 P, 6 D, 6 F, 
e G odd-parity LS'-terms, the wavefunction expansion consists of 1102 CSFs. 

The wavefunction expansion for the lower state in the case of the 2s — 5/ transition 
array is constructed in a similar way by replacing all the configurations containing 4/ 
(/ = 0, 1,2,3) electrons with configurations containing 5/ electrons. Furthermore, a complex 
of configurations with one electron in the 4/ shell is added. This extends CI expansion 
by including odd-parity configurations 2p 6 3/ 2 4/' (/ = 0, 1, 2; /' = 1, 3) and 2p 6 3l3l'Al" 
(/ = 0, 1; /' = 1, 2; /" = 0, 1, 2, 3). The same (as for the n = 4 lines) LS'-terms are considered, 
and the CI wavefunction expansion gives rise to 716 and 1193 CSFs for the J = 1/2 and 
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J = 3/2 levels, respectively. 

While considering the upper state of Fe 13+ , the CI wavefunction expansion for 2s — 3/, 4/ 
lines is constructed of several configuration complexes, representing configurations with 
one or two electrons virtually excited from the ground configuration. We have included 
configurations 2p 6 3s 2 3d, 2p 6 3s3p 2 , 2p & 3s3d 2 , 2p 6 3p 2 3d, 2p 6 3d 3 with three outer electrons 
in n = 3 shell, the even-parity configurations 2p 6 3l 2 Al' (Z = 0,1,2;/' = 0,2), 2p 6 3Z3Z'4Z" 
(I = 0, 1; V = 1, 2; I" = 0, 1, 2, 3), 2p 6 3lAl' 2 (Z = 0, 2; V = 0, 1, 2, 3) with one electron in outer 
n = 4 shell, 2p % 3lAl'Al" (Z = 0, 1, 2; V = 0, 1, 2; I" = 1, 2, 3) with two electrons in n = 4 shell. 
On top of that, we have extended the CI expansion by adding the even-parity configurations 
with a vacancy in the 2s shell: 2s2p 6 3l 2 3l' 2 (I = 0,1; Z' = 1,2), 2s2p 6 3s3p 2 3rf, 2s2p 6 3s3d 3 , 
2s2p 6 3p 4 , 2s2p 6 3d 4 , 2s2p 6 3l 2 3l'Al" (I = 0,1; V = 1,2; I" = 0,1,2,3), 2s2p 6 3/3/ ,2 4/" 
(I = 0, 1; V = 1, 2; I" = 0, 1, 2, 3), 2s2p 6 3/ 3 4Z / (Z = 1, 2; Z' = 0, 1, 2, 3), 2s2p 6 2s2p3dAl (1 = 1,3), 
2s2p & 3l 2 Al' 2 [l = 0,1, 2;V = 0,1,2,3), 2s2p 6 3l 2 4l'Al" (Z = 0, 1,2;Z' = 0,1; I" = 2,3), 
2s2p 6 3s3rf4Z 2 (Z = 0,1,2,3), 2s2p 6 3l3l'Al"Al'" (Z = 0,1; I' = 1,2; I" = 0,1,2;Z W = 1,2,3). 
We consider the same LS'-terms as in the transitions from the 2p shell. Consequently, the 
numbers of CSFs included in the CI wavefunction expansion are 1695, 2865 and 3270 for 
the J = 1/2, 3/2, 5/2 levels, respectively. 

While constructing the upper-state CI wavefunction expansion configuration set for 
the 2s — 5Z transition array, we apply the same principles as in the n = 4 case but replace 
the configuration sets containing 4Z electrons with those containing 5Z electrons. That 
expansion is extended with a set of even-parity configurations with a valence n = 4 electron 
2p 6 3Z 2 4Z' (Z = 0,1,2;Z' = 0,2) and 2p 6 3l3l'4l" (Z = 0,1; I' = 1,2; I" = 0,1,2,3). This set 
generates 1751, 2954 and 3364 CSFs in the CI wavefunction expansion for the J = 1/2, 3/2 
and 5/2 levels, respectively. 
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2.3.3. Transition data 

Using the above methods to generate the CI wavefunction expansions, we obtained 
wavefunctions to produce photoexcitation line data for the Fe XIV ion. Calculated line data 
cover the wavelength region A = 10.571 — 14.432 A for 2s — 31 transitions (a total of 617 
lines), A = 9.298 - 11.707 A for 2s -4/ transitions (3194 lines) and A = 8.890 - 10.861 A for 
2s — 51 transitions (147 lines). For transitions from the 2p shell, the wavelength range was 
A = 12.341 -16.626 A for excitation from the n = 3 shell (1343 lines), A = 10.141 - 13.642 A 
for the n = 4 shell (7425 lines), and A = 9.685 - 12.438 A for the n = 5 shell (8640 lines). 

Table [3] lists the atomic data for the Fe xiv UTA lines. Transition wavelengths A and 
absorption oscillator strengths fij for excitation from the 2s and 2p shells to the outer 
n = 3, 4, 5 shells are given for lines with / > 0.1. A more extensive data set is presented in 
the online version of this table, where we provide results for all lines with / > 0.0001 (see 



Sect. 3.2 for more details about data trimming) 



In Fig. [3] we compare the strongest lines from our calculation with the results derived 



using the MBPT approach of Gu et al. (2006) and with the HULL AC relativistic data from 



Behar et al. (2001). There is reasonable wavelength agreement for all three sets, as found 



for Fexvi and Fexv. The only exception is the line with A = 16.033A from Behar et al. 



(2001), with the larger values of oscillator strengths ^2 fij found by these authors explained 



in Sect. 2.1.3 A wavelength disagreement between our results and the MBPT data was 



also discussed in that section. 
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Table 3. Transitions 2p — rijd and 2s — rijp from the L shell to the rij = 3, 4, 5 shells in 
Fe 13+ . Column n^i denotes the initial shell, rij the final shell for optical electron transition, 
Qi and Qj the statistical weights for the initial and final levels, A the line wavelength (in A) 

and fij the absorption oscillator strength. 
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14 15 
Wavelength [A] 



Fig. 3. — A comparison of the wavelengths A and absorption oscillator strengths /j,- for the 



strongest lines in Fexiv. Dots represent our data, crosses are the results of Gu et al. (2006) 



and stars those from Behar et al. (2001). The data at A = 15. 5A represent the 2p — 3d 



UTA, those at A = 14. 3A represent the 2s — 3p UTA, and those at A = 12. 9A represent the 



2p — Ad UTA. The Behar et al. (2001) data are given for an /-value averaged wavelength 
A av and the total oscillator strength ^ f%j of the line group. See the electronic edition of the 
Journal for a color version of this figure. 
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3. Incorporation into Cloudy 
3.1. Data sources 

The UTA data sets now used in Cloudy are summarized in Table |4| The Opacity 



Project (OP; Badnell et al. 2005) produced complete data, for the K shell as well as the 
LI and L2 shells where appropriate, for ions with twelve or fewer bound electrons. These 
include the auto-ionization branching ratio which we consider in calculating the ionization 



rate. The OP data are marked "B" in this table. Gu et al. (2006) provide all stages of 
ionization for Fe but includes only excitations from the 2p shell, denoted by "G" in the 
Table. These data sets have been included in Cloudy since shortly after their original 
publication dates. 

The data presented in this paper are now being incorporated into Cloudy and are 
denoted by "K" in Table @ Two ions, Fe 14+ and Fe 15+ , are present in both the current 
data and the OP. We use the OP data for excitations from the K shell and the current 
results for L shell excitations. Our new data do not include autoionization rates, so these 



were copied from Gu et al. (2006) for the 5 most important lines of FeXV and FeXVI. For 
other lines we assume that the branching ratio for ionization is 100%. In the default setup 
of the code, we only use transitions from the ground level for all data sources. See also the 
discussion in 5 14.11 



It is clear from Table [4] that many ions simply have no current data. For example, 
there are no K shell data for Fe 13+ . These are all major requirements for accurate spectral 
simulations. 



- 26 - 



00 



> — < 



O 



to 
10 



a 



pq 

a a 



m ea a pq 



K X 



pq 

pq pq 
pq pq pq 
pq pq pq pq 
pq pq pq pq pq 

CD 



pq 
pq 
pq 
pq 
pq 
pq 



pq 
pq 
pq 
pq 
pq 
pq 
pq 
pq 



pq pq O Z O fa £ 



pq 
pq 
pq 
pq 
pq 
pq 
a 
pq 
pq 



pq pq pq pq 
pq pq pq pq 
pq pq pq pq 









pq 


a 






pq 


pq 


pq 




pq 


pq 


pq 


a 


pq 


pq 


pq 


pq 


pq 


pq 


pq 


pq 


pq 


a 


pq 


pq 


pq 


pq 


a 


pq 


pq 


pq 


pq 


a 


pq 


pq 


pq 


pq 


a 


pq 


pq 


pq 


pq 


a 


pq 


pq 


pq 


pq 




pq 


pq 


pq 






pq 


pq 








pq 











a 

a a 

a a a 

a a a 

a a a 

a a a 

a a a 

a a a 

a a a 

a a a 
a a 
a 



a a a a a 

a a a a 

a a a 

a a 



a 



a cq o < W 



C3 M 



H > O 



-27- 



o 

a 
+3 
pi 
o 
O 



03 

03 03 

cq m m 

cq m m 

m m m 

eq M pq 

m m m 

to n to 

03 m 03 



03 

03 03 

03 03 03 

03 03 03 

03 03 03 

03 03 03 

03 03 03 

03 03 03 

03 03 03 

03 03 03 
03 03 
03 



03 



03 



03 



03 



03 
03 « 

o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 



o 

OS 



3 a 

•Z, O N 



-28- 



3.2. Data trimming 

The newly computed data sets are quite extensive and include a large number of 
exceptionally weak lines. Figure [4] shows the absorption line oscillator strength fi u plotted 
against the normalized sum of all lines with / > fi >u . Careful examination of such data 
shows that the summed oscillator strength, which affects the UTA ionization rate, has 
converged to within 0.57%, 0.15%, and 0.08% of the total for Fe 13+ , Fe 14+ , and Fe 15+ , when 
all lines with / > 10~ 4 are included. Weaker lines have little effect upon the ionization 
rate because the line-center opacity is smaller than the continuous photoionization opacity. 
Accordingly, we only consider lines with / > 10 -4 , although we retain all data for future 
flexibility. For the case of Fe 13+ this reduces the number of lines from 80894 to 754, for 
Fe 14+ from 2738 to 103, and for Fe 15+ from 595 to 137. 



3.3. Damping parameters and line broadening 

The Voigt function, which describes the line profile including both thermal and natural 
broadening, is given by 

a [ + °° exp(- j/ 2 ) 

H{a,x) = - dy (3) 

n J-oo {x-y) + a 2 

where x is the displacement from line center, measured in terms of the Doppler width 

The damping parameter a is the ratio of natural (radiation damped) to thermal line widths 

4nAu Dop 

and the natural broadening width is given by the damping constant 7 which is expressed as 
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Fig. 4. — Absorption line oscillator strength fi :U plotted against the sum of all lines with 
/ > fi jU . The y-axis has been scaled by the sum of the oscillator strengths. The integrated 
absorption is well converged when only lines with / > 10 -4 are included. 

a sum of radiative A r and autoionization A a rates: 

7 = E^ + E^' ( 6 ) 

Ku l'<u 

H(a, x) is normalized so that its integral over x is y/n. 

The Voigt function must be evaluated to account for line self-shielding when finite 
column densities are encountered, while the total radiative decay ratio out of the upper level 
is needed to derive the damping constant 7. This lifetime is the sum of the autoionization 
and radiative decay rates. Our new data do not include calculations of auto ionization 
rates (except for the 5 lines mentioned earlier). However, in order to calculate the damping 
constant 7 we need to assume a value, so we adopted the value, 10 15 s _1 , taken from Opacity 
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Project results for Fe XV and Fe xvi. The other data sources include auto ionization rates, 
which we adopt. 

Figure [5] shows the resulting damping parameters, plotted as a function of transition 
wavelength, for a gas kinetic temperature of 10 4 K. UTA transitions are strongly damped, 
often with a> 1, because of the rapid autoionization rate. For comparison, the damping 
parameter for a strong UV line, such as H I La, is a ~ 4 x 10 -4 . It should be noted that the 
high damping parameter a is a consequence of the inner-shell nature of these excitations 
regardless of their being in UTAs. 

While simple approximations to the Voigt function may be sufficient for UV - optical 
spectroscopy, with small a, the large range in damping parameter that occurs when UTAs 
are considered make it important to use implementations of the Voigt function which 



are accurate for all a. We adopt the routine provided by Wells (1999), combined with a 



specially designed faster routine for a < 0.1, which provides results with a relative accuracy 
of 1 in 10 4 , and has been confirmed to pass all the test cases given by Zaghloul & Ali (2011 ). 



Stellar atmosphere texts (Rutten 2003; Mihalas 1978) often focus on results derived 



for approximations to the Voigt function which, while valid in the original context, are not 
accurate for the full range of frequency and a needed to include UTA transitions. Figure [6] 
shows the function for a typical strong UV line such as H I Lya (a ~ 10~ 4 ) and a strongly 
damped UTA transition (a ~ 10 3 ). We see that at line center H (a,x) ~ (1 + a) -1 , and 
that the core of the line is roughly x ~ (1 + a) wide. 

There are a large number of lines which strongly overlap to produce the observed 
UTA features. The effects of mutual line shielding due to overlap must be included if the 
ionization rate is to be properly computed. Line overlap is treated using a combination of a 
coarse, low resolution continuum, and a fine continuum, with resolution sufficient to resolve 



lines, as described in Shaw et al (2005). 
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Fig. 5. — Damping parameters as a function of the line wavelength for our complete set of 
UTA lines. 




Fig. 6. — The Voigt function H(a,x) for two values of the damping parameter a. 
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4. Application to AGN 

4.1. Fe XIV as a density indicator 

Our Fe XIV data include both fine-structure levels within the ground term, namely 
2 Py 2 an d 2 -^3/2 • These are separated by 2.33741 eV and produce the famous "green" Fexiv 
coronal line at 5303 A. There are approximately the same number of UTA lines originating 
from either level in the overall configuration. 

The critical density of the excited J = 3/2 level is ~ 3 x 10 9 cm -3 at 10 5 K. The 
temperature T = 10 5 K is approximately that of the formation of the Fe ions discussed 
in current paper. At considerably lower densities most of the population will be in the 
J = 1/2 level, while for high temperatures and densities the levels will be populated 
according to their statistical weight and most will be in the excited level. If the absorption 
characteristics are different in these two different limits then the UTA lines could be used 
as a density indicator. 

Figure [7] (right frame) shows how the profile of the A15.5 A Fe XIV UTA changes with 
density. The calculation assumed a total hydrogen column density of 10 21 cm -2 , a solar Fe 
abundance, and two hydrogen densities, n-^ = 10 5 cm -3 , sufficiently low for all populations 
to lie in the lower J = 1/2 level, and nn = 10 15 cm -3 , high enough for the levels to be 
populated according to their statistical weights. Significant differences are present. The 
features originating from the configurations which include the excited J = 3/2 level are 
not present in the low density profile. At the high density limit, the ratio of the excited 
and lower level population is 2:1. Therefore, new features connected to the J = 3/2 level 
(e.g. A15.35A, A15.40A, A15.65A) become evident, whereas features originating from the 
ground J = 1/2 level (e.g. A15.30A, A15.43A, A15.47A) become less pronounced. The 
total absorption does not change because the sum of the absorption oscillator strengths in 
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Fig. 7. — UTA spectrum of Fe XIV transmitted through a cloud with a column density of 
N(H)= 10 21 cm -2 and a solar Fe abundance. The effects of different densities upon the shape 
of the Fe XIV UTA are shown, with the low density case in blue having nearly all populations 
in the lower J = 1/2 level, while the high density case has level populations within the ground 
term proportional to their statistical weight. This calculation only included Fe XIV. 

this wavelength region is approximately the same both for J = 3/2 and J = 1/2 levels, 
although different transitions are produced by each level. 

Figure [7] (left frame) shows the effects of density in the area of the A13.0A UTA lines 
in Fe XIV. It demonstrates very similar behavior to the case of the A15.5A UTA shown in 
the right panel, with J = 3/2 lines appearing at higher densities and J = 1/2 lines being 
prominent at low-density limit. 

Future generations of X-ray spectrometers may be able to use such differences to 
measure the density of the absorbing gas. This would make it possible to deduce the 
location of the gas. Nevertheless, improved accuracy of the theoretical wavelengths would 
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be needed for such a test to be definitive. 

Absorption line spectroscopy generally has a degeneracy introduced by the fact that the 
ionization distribution determines an ionization parameter, the ratio of the flux of ionizing 
photons to gas density. Literature on this type of analysis is extensive, and summarized 



by Chakravorty et al. (2009). The ionization itself is not sensitive to either flux or density. 
By measuring the density and ionization parameter, the flux could be deduced, which then 
leads to the determination of the source - cloud separation. 

Tests show that the density effects do not have a significant effect on the ionization 
of the gas because the total absorption oscillator strengths are similar. In the current 
implementation we assume the low density limit in computing the effects of UTA transitions 
and the resulting spectra. 

The ground terms of Fe 14+ and Fe 15+ are simpler, having only one level, so such 
distinctions do not occur. 



4.2. Probing the thermal stability of the warm absorber 



The nature of the absorbing and emitting clouds in AGN is a long-standing problem 



(Osterbrock & Ferland 2006 Chakravorty et al. 2009). By analogy with the local ISM 



several gas phases are thought to be present at one location, each in equilibrium with the 
radiation field, but having different levels of ionization and a range of temperatures, but the 



same gas pressure. This has been discussed extensively in the literature (Krolik et al. 1981 



Hess et al. 1997; Krolik & Kriss 2001 Komossa & Mathur 2001 Reynolds & Fabian 1995). 



The ionization, temperature and spectrum of clouds that can exist is determined by 
the type of stability analysis shown in Figure |8j This shows the familiar "S curve" , which 
is computed for an optically thin cell of gas exposed to the AGN radiation field. Thermally 



- 35 - 



stable phases have positive slope, while regions with negative slope are unstable so that gas 
will only exist in these regions for a short time. The shape of the S curve determines the 
properties of clouds which are thermally stable and so live long enough to contribute to the 
observed spectrum. 

10 7 i j 
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10 5 



-0.3 -0.2 -0.1 

log[U/Tx 10 4 [K-']] 



0.1 



0.2 



Fig. 8. — A newly-computed AGN thermal stability curve. The default AGN SED described 
by Chakravorty et al. (2009) and a hydrogen density of 10 8 cm -3 are assumed. Thermally 



stable regions are shown by the thicker lines. 



Recent S-curve calculations have focused on how properties of the AGN might affect its 
shape, and so determine which clouds might exist. The literature on this topic is vast, and 



recent examples include Holczer, Behar, & Kaspi (2007), Holczer, Behar, & Arav (2010) 



and Holczer & Behar (2012), who examine observational determinations of the stable points 



on the S curve, Chakravorty et al. (2008) who show the effects of updated dielectronic 



recombination rate coefficients, Chakravorty et al. (2009), who do a systematic stability 
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analysis and and Chakravorty et al. (2012), who discuss the influence of the SED on the 
cloud stability. 

Here we show the effects of the improvements discussed in this paper, using version 



C13 of Cloudy dFerland et al. 2013). Some other recent improvements to the atomic data, 



which focus on the atomic models used to compute the cooling, are described by Lykins et 



al. (2013). Figure H shows the thermal stability curve resulting from the improved atomic 



physics in this paper. Many model parameters are taken from Chakravorty et al. (2009). 



The x-axis is the ratio of the dimensionless ionization parameter U, defined as the ratio of 
ionizing photon to hydrogen densities, to the gas kinetic temperature, while the latter is the 
y-axis. Thermally stable regions, those with positive slope, are shown as the heavier lines. 



In a series of papers Holczer, Behar, & Kaspi (2007), Holczer, Behar, & Arav (2010) 



and Holczer & Behar (2012) infer from observations of column densities of certain ions that 



temperatures between 4.5 < log T < 5.0 are missing and ascribe it to thermal instability. 
On the other hand, this region is stable in the current calculation with its assumed 
parameters. The shape of the stability curve is affected by several other ingredients besides 
the atomic data. The composition assumed, and the form of the SED, also change it. This 
suggests that the presence or absence of stable gas could be used to infer the SED or gas 
composition, among other properties of the AGN. 

Figure [9] shows the distribution of ionization stages of Fe as a function of the ionization 
parameter U. The ions discussed in this paper peak in the range log U ~ 1 — 2. These 
probe the upper bound of the low-T branch of thermally stable gas, the lowest unstable 
region, and the low-T part of the middle stable branch. Horizontal lines in the upper part 
of Figure [9] indicate the regions where the gas is stable. 

The atomic data presented in this paper affect the details of the transition between the 



cool and middle stable branches. Figure 10 shows the ratio of UTA to total ionization for 
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Fig. 9. — The Fe ionization distribution across Figure [8j where the gas kinetic temperature 
is indicated by the solid line increasing from left to right. Horizontal lines at the top of the 
figure mark the thermally stable regions of Figure |8j The Fe ions considered in this paper 
probe the high temperature end of the cool branch in Figure [8j unstable regions, and the 
low temperature end of the middle stable branch. 

several charge states of Fe, where the filled circles include the data presented in this paper, 
while the crosses represent the older data. The larger number of lines derived here results 
in a larger UTA ionization rate, as shown in the Figure. By UTA ionization we mean the 
effect of autoionization following inner-shell photo-excitation. UTAs have the greatest effect 
on the ions which occur around the transition from the cool to middle stable branch, and 
their physics affects the details of this transition. 



Finally, in Figure 11 we show the spectrum produced by an intervening cloud near the 
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Fig. 10. — The ratio of the UTA ionization to total ionization rate for several Fe charge 
stages. Filled circles include the results presented in this paper, while the crosses represent 
the older data. These are for the conditions occurring across Figure [8] Rates are evaluated at 
the illuminated face so line self-shielding is not important. UTA ionization is the dominant 
ionization process for the ions considered in this paper. 

upper range of the cool stable branch. This cloud has solar abundances, a column density 
of N(H) = 10 21 cm" 2 , and an ionization parameter of log U = 1.25. A portion of the coarse 
continuum which is used for continuum radiative transfer and atomic rates is shown in the 
upper panel. The assumed SED, which includes the "Big Bump" from the central accretion 
disk and a non-thermal X-ray power law, has many absorption features superimposed. 
Emission lines are also produced by the cloud but are weak in this portion of the spectrum 
due to the low gas temperature. 

The lower panel shows a small portion of the fine continuum in the neighborhood of the 
~ 15 A UTA feature. Two curves are shown, with the solid line using the results presented 
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in this paper, and the dashed line employing the other data sources summarized above. 
Significant differences are present. 
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Fig. 11. — The spectrum of a cloud lying along the line-of-sight to an AGN. Cloud parameters 
are indicated in the upper panel, which shows the transmitted coarse continuum. The lower 
panel shows a small portion of the fine continuum with predictions using the atomic data 
summarized in this paper in black, and with older UTA line data in green. 
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5. Conclusions 

This paper has summarized advances in the atomic data needed to simulate conditions 
in a photoionized plasma. All of these improvements are now in the development version of 
the plasma simulation code Cloudy, and will be part of the next major release. The specific 
advances are the following: 

• We have calculated a large set of atomic data for UTA lines using the CI method 
implemented in the Civ3 code. This new dataset substantially extends previous 



atomic data of Kisielius et al. (2003) in two ways. First, the El transitions from the 
inner 2s shell are determined. Second, the data sets for Fe XIV, Fe XV and Fe XVI 
include lines 21 — 41, 5/ in addition to the earlier determined lines 2p — 3/. 

We have incorporated this large set of new UTA data into the spectral simulation 
code Cloudy, and applied it to problems in AGN. These improvements will be part of 
the next major release of the code. 

We summarize our data sources for UTA transitions. There are still missing data, even 
for very important ions. These should be a priority for new theoretical calculations. 

We summarize how our data, which were computationally very expensive to undertake, 
compare to simpler calculations. Line wavelengths differ (although insignificantly) 
due to differences in the computed level energies, but the transition rates are in good 
agreement. 

The UTA lines are often strongly damped, many having damping parameters a ^> 1. 
We have improved the form of the Voigt function used by Cloudy to handle such 
strongly damped lines. 

We show how the Fe XIV UTA at A15.5A can be used to measure the density of the 
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gas, or identify whether the density is significantly above or below 10 9 cm 3 . Such 
measurements would help determine the location of the warm absorber in AGN. 

• The total ionization rate is increased by roughly 30% with the new set of UTA data, 
which have far more lines. This changes the ionization of the gas and alters the 
thermal properties of photoionized gas exposed to a typical AGN SED. 

• We present a newly computed thermal stability "S curve" using the new data. We 
show that the Fe ions considered in this paper are produced in the warmer parts 
of the cool thermally stable branch, an unstable region, and in cooler parts of the 
middle stable branch. As a result, these lines probe the portions of the S curve which 
determine which cloud parameters can persist. Future work will investigate the effects 
changes in these regions have upon the observed spectrum. 
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